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I. INTRODUCTION 



It has been known that two neutral atoms or molecules separated by a distance which is rather small but much larger 
than the atomic dimensions interact through the van der Waals force. As was shown by London, the van der Waals 
force arises in second order perturbation theory from the dipole-dipole interaction. It is caused by the dispersions of 
dipole operators, i.e., by quantum fluctuations Being a quantum phenomenon, it depends on the Planck constant 
ft. In consequence of interatomic interactions, the van der Waals force acts also between an atom and a macroscopic 
body and between two closely spaced macroscopic bodies. If the separation between two atoms, an atom and a 
macroscopic body or between two macrobodies is sufficiently large, so that the retardation of the electromagnetic 
fluctuating interaction contributes significantly, the van der Waals force depends on both ft and the velocity of light 
c. In this regime the force is usually labeled as Casimir Q or (in the case of atom-atom or atom-wall interaction) 
Casimir-Polder Q force. Van der Waals and Casimir forces play an important role in the interaction of a single layer of 
graphite (hereafter, graphene) with macroscopic bodies (a material plate or a semispace) and with microparticles (an 
atom or a molecule). The physics of interactions between graphene and carbon nanostructures Q with macroscopic 
bodies and microparticles is significant for the understanding of layered systems, bundles of nanotubes or metallic 
nanowires and absorption phenomena. Special attention was attracted to this problem after the proposal of Ref. [f| 
to use the single-wall carbon nanotubes for the purposes of hydrogen storage. Subsequently controversial results on 
this subject were obtained Q. Especially, the microscopic mechanisms underlying the absorption phenomenon remain 
unclear. 

Theoretically the interaction of hydrogen atoms with graphite sheets and carbon nanotubes has been studied 
in Ref. 0] using the density functional theory and in Refs. [1, Q using the nonrelativistic perturbation theory for 
degenerate levels of a two-level atomic system. The adsorption of fullerene molecules on graphite was considered 
in Ref. [Icl ] demonstrating that the interaction is basically proportional to the number of neighboring atoms. Most 
of theoretical work on the van der Waals force in layered structures and between carbon nanotubes was done using 
the density functional theory [ll|, El, E!> EH EM EH • It is known, however, that in some cases density functional 
theory (especially when used with linear-density approximation) gives not enough precise description of van der Waals 
interactions [l7[- In Ref. [H| the Lifshitz theory of the van der Waals and Casimir force was extended to the case 
of a microparticle interacting with a plane surface of an uniaxial crystal and with a multi-wall carbon nanotube. It 
was shown that the position of a hydrogen atom inside the multi-wall carbon nanotube is energetically preferable as 
compared to the outside position. 

It is common knowledge that the Lifshitz theory [h| [2(| provides the fundamental description of the van der Waals 
and Casimir interaction between two macroscopic bodies and between a microparticle and a macroscopic body. In the 
framework of this theory the interaction energy and the force are expressed in terms of the dielectric permittivity of 
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macrobodies or the dynamic polarizability of a microparticle. Recently, the Lifshitz theory was successfully applied to 
the interpretation of precision measurements of the Casimir force 2^. 12, US [H]j25| and to the calculation of atom- 
wall interaction in connection with Bose-Einstein condensation [26|, [27| . In Ref. [18| the Lifshitz theory was adapted 
for the case of sufficiently thick multi-wall carbon nanotubes using the idealized description of the wall material by 
dielectric permittivity. This idealization is, however, not applicable for the description of single-wall nanotubes which 
narrows the applicability of the standard Lifshitz theory. 

In the present paper we use the description of graphene in terms of the two dimensional free electron gas [28j 
in order to extend the Lifshitz theory of the van der Waals and Casimir interaction to the case of carbon systems. 
Graphene is considered as an infinitesimally thin positively charged flat sheet, carrying a continuous fluid with some 
mass and negative charge densities. This sheet is characterized by some typical wave number £1 determined by the 
parameters of the hexagonal structure of graphite. In Refs. [1^, [3(| the interaction of the electromagnetic oscillations 
with such sheet was considered and the normal modes and reflection coefficients were found. The van der Waals and 
Casimir interaction between the two parallel plasma sheets was described in Ref. [3l| using a Lifshitz-type formula. 
The important distinctive feature of this formula is that it does not use the concept of the dielectric permittivity of 
the sheet. 

Here we obtain Lifshitz-type formulas for the van der Waals and Casimir interaction of graphene and material 
semispace or a plate. In so doing the semispace or the plate are described using the dielectric permittivity and 
graphene is considered as a plasma sheet. The numerical computations of the interaction energy and the force acting 
between graphene and thick plates made of different materials arc performed. The Lifshitz-type formula for the 
interaction of graphene with a microparticle (an atom or a molecule) is also obtained. This formula is applied to 
the computations of the interaction of graphene with hydrogen atoms and molecules. The approximate formulas for 
the interaction of a single-wall carbon nanotube with a material plate are presented when the nanotube is in close 
proximity or far away from the plate. 

The paper is organized as follows. In Sec. II the Lifshitz-type formulas are obtained for the interaction of graphene 
with a semispace or a plate. Sec. Ill contains the results of numerical computations using these formulas for metal 
and dielectric plates. In Sec. IV the interaction of graphene with a hydrogen atom or a molecule is considered. Sec. V 
is devoted to the interaction of a single-wall carbon nanotube with a metal or dielectric plate. Sec. VI contains our 
conclusions and discussion. 



II. INTERACTION OF A GRAPHENE WITH A SEMISPACE OR A PLATE 



We consider the van der Waals and Casimir interaction of a graphene occupying the xy-plane, z = 0, with a material 
semispace or a plate of thickness d. The separation distance between the boundary plane of the semispace (plate) and 
graphene is a. The dispersion interaction of the two plane parallel bodies (plates or semispaces) labeled by the upper 
indices 1 and 2 with the electromagnetic oscillations can be described in terms of the reflection coefficients TE 

(2) 

and Trpj^ TE for two independent polarizations of electromagnetic field (transverse magnetic and transverse electric). 
In so doing the van der Waals and Casimir interaction energy per unit area is given by the Lifshitz formula 0, [2(| Hl| 



In 1 1 - r (1) r (2) e^ 2aq 



(1) 



Here k±_ is the magnitude of the wave vector component perpendicular to the z-axis (i.e., lying in the plane of the 
plates), £ is the frequency variable along the imaginary axis (u> = it;), and 



£ 2 

(2) 



r-2 



Equation (TlJ is applicable at not very high temperatures (not too large separations). As an example, at room 
temperature (T = 300 K) thermal corrections to Eq. (fTJ) are negligible up to separations of about 1 fim. If the 
temperature is higher (separation is larger) Eq. (fT]) is simply generalized by changing the integration with respect to £ 
for the summation over the discrete Matsubara frequencies. The general derivation of Eq. (fT]) for arbitrary reflection 
coefficients can be found in Ref. (2lj . 

The commonly accepted approach (see, e.g., 0, [2(J HH, [H, [24], [2^, [26|, [2?| ) uses the Fresnel reflection coefficients 
of the semispaces and plates expressed in terms of the frequency-dependent dielectric permittivities. Let the semispace 
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labeled 2 be made of isotropic material and be described by the dielectric permittivity e(u). In this case the reflection 
coefficients are [H, M, EH [H H [H [H 

_(2) =r (2) f , fe ,^_ g(*Qg~fe 

r TE = r TE,s(£' ^J-) = ^ 

where 



If the second body is a plate of finite thickness d, Eq. ([3]) should be replaced by [H| 

(2) = (2) i \ £ 2 (»C)g 2 ~ fc 2 

r ™ ~ r ™^ ? ' 1J e 2 (*0<7 2 + ^ 2 + 2gte(i0 coth(fcd) ' 
(2) (2) , k 2 -g 2 

" TE,pVS j fc i J - fc2 + ?2 + 2gfc CQth(fcd) " W 

If the first body is also a semispace or a plate, its reflection coefficients are obtained from Eqs. and ^ by the 
replacement of the upper indices 2 for 1. 

In our case the first body is graphene and it cannot be described macroscopically it terms of the dielectric per- 
mittivity. The van der Waals and Casimir interaction of a semispace or a plate made of isotropic material with the 
graphite plate containing a few hexagonal layers was considered in Ref. [181 ] using the reflection coefficients for a 
uniaxial crystal. In the case of a graphite plate of thickness d they take the form 



' TM — r TM,p 



M)zM)<l 2 + + 2qk zy /e x {i£)eM) coth(fc*d) ' 
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r w= r w (tk,)- x ~ q a) 

where £x(u)) — £ y (co) and £ z (oj) are the graphite dielectric permittivities in the directions x,y and z, respectively, and 



/•' , = \/ k\ + £ X % , k z = Jki + £ z (iO % . (6) 



The number of layers in the graphite plate should be sufficiently large that the macroscopic description in terms of 
the dielectric permittivities is applicable. 

The reflection coefficients for a sheet of graphene cannot be obtained from Eq. ([5]) in the limit d — > [in fact, the 
coefficients (J5J) go to zero when d vanishes]. The reason is that the case of "thin" plate implies that d/a is sufficiently 
small, whereas d should be large enough for the validity of the macroscopic description in terms of e. In Sec. Ill the 
correlation between the macroscopic description of a graphite plate by means of Eq. ([5]) and the case of graphene will 
be clarified. 

A single plane hexagonal layer of graphite, i.e., graphene, can be described as infinitely thin plasma sheet where 
the 7r-electrons are treated as a continuously charged fluid moving in an immobile, overall neutralizing back grou nd 
of positive charge. These plasma sheets have been considered in Ref. [28[ and, more recently, by Barton (29l. |30| in 
connection with the Casimir effect for a fullerene and for a single base plane from graphite. In Ref. 31] they were 
used to calculate the van der Waals and Casimir interaction between the two parallel graphenes and a Lifshitz-type 
formula was obtained for their interaction energy. 

The plasma sheet model assumes the 7r-electrons of the carbon atoms to be described by a negatively charged fluid 
confined to a plane and having the two dimensional displacement vector R(x, y) exp(— iuit), where we introduced at 
once the usual harmonic time dependence. (In [29|, [3(| this vector is denoted by £ and the wave number of the sheet 
by q.) The fluid has a surface charge density ne and a surface mass density nm where e and m are the electron charge 
and mass, respectively. For the hexagonal structure of carbon layers there is one 7r-electron per atom [32| resulting in 
two 7r-electrons per one hexagonal cell. This leads to 
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where I = 1.421 A is the side length of a hexagon. 

The fluid provides a source for the Maxwell equations with surface charge and surface current densities, 

a = — neVt • -R, j = —iuineR, (8) 

where the operator V< acts in the tangential direction to the sheet [here and below a, j and fields E and B depend 
on coordinates; their dependence on time is obtained through the multiplication by the common factor exp(— itot)]. 
The Maxwell equations read 

V • E = AnaS(z), V x E - — B = 0, (9) 

c 

V-B = 0, V x B + — E = — jS(z). 

c c 

By the integration of Maxwell equations across the sheet, we obtain the matching conditions on the tangential and 
normal components of the fields [3(| > 

c 2 

E t .2 - E t ,i = 0, E z , 2 - E z<1 = 2fi-jV t • E t , (10) 
B za - B zS = 0, B t 2 - B t ,i = -2Xl-z x E t . 

UJ 

Here z — (0, 0, 1) is the unit vector pointing in z-direction, 

77 P 2 

fl = 2vr = 6.75 x 10 5 m" 1 (11) 

mc 2 

and 77 is defined in Eq. ((7|). The quantity f2 is the main characteristic of graphene in the model under consideration. 
The value in Eq. pT|) corresponds to the frequency = c51 = 2.02 x 10 14 rad/s (as compared, for instance, to the 
plasma frequency of gold uj p — 1.37 x 10 16 rad/s). 

We remark that now the matching conditions (JTOJ) together with the Maxwell equations © outside the surface, 
i.e., without the delta functions on the right-hand sides, provide the complete description of the interaction of the 
electromagnetic field with the plasma sheet. This implies, for instance, that all components of the field satisfy the 
usual Poisson equations, 

A + ^E = Q, (A + ^B = 0. (12) 

After the separation of variables x and y in Eqs. (jlpp and (|12[) we arrive at a one-dimensional scattering problem in 
the z-direction [21[ . The solution of this problem leads to the following reflection coefficients on the graphene plasma 
sheet taken at the imaginary frequency axis [30| . 

4 1 ^r« g (C,fex) = ^. (13) 

Equations (fT]) and ©, (TIB"]) or, alternatively, (fTJ) and ([3]), (TIB")) allow to calculate the energy of the van der Waals 
and Casimir interaction between a graphene and a semispace or a plate made of some usual material described by 
the dielectric permittivity e(u>). From Eq. (TT]) it is easy to obtain the Lifshitz-type formula for the van der Waals and 
Casimir force per unit area acting between graphene and material semispace or a plate, 



r (l) r (2) r (l) r (2) 

'tm'tm i 'te'te 



olaq _ J 1 ) r ( 2 ) „2aq _ r W r ( 2 ) 
'TM'TM e 'TE'TE/ 



(14) 



Equations (fT]) and (TT4]) can be used to obtain the interaction energy and the force between an atom (molecule) and 
a graphene sheet (see Sec. IV) and between a carbon nanotube and a material wall (Sec. V). In the next section these 
equations are applied to compute the interaction of graphene with metal and dielectric walls made of real materials. 



5 



III. COMPUTATIONS OF THE VAN DER WAALS AND CASIMIR INTERACTION BETWEEN 
GRAPHENE AND METAL OR DIELECTRIC PLATE 



Here we apply Eqs. (Q]) and (fT4|) to calculate the van der Waals and Casimir interaction energy and force acting 
between graphene and thick plates (semispaces) made of different materials. The reflection coefficients TE for 

graphene are given in Eq. (fT"3|) , and the reflection coefficients TE for a semispace are presented in Eq. ([3]) . All 
parameters in Eq. (|13[) are specified. As to Eq. ([3]), it depends on the dielectric permittivity along the imaginary 
frequency axis. In our computations we consider metal and dielectric semispaces made of Au and high resistivity Si, 
respectively. High precision results for the dielectric permittivities of these materials along the imaginary frequency 
axis were obtained in Ref. [33[ by means of the Kramers-Kronig relation 

* Jo u + £ 

The imaginary part of the dielectric permittivities at real frequencies was taken from the tabulated optical data [34j . 

The computational results for the van der Waals and Casimir energy density E(a) normalized to the Casimir energy 
density in the configuration of the ideal metal plates, 

„ , ^ tt 2 he , . 

*«— TO?' < 16 > 

are shown in Fig. 1(a) as a function of separation. The solid and dashed lines are related to the interaction of graphene 
with Au and Si, respectively. In Fig. 1(b) the analogous results for the van der Waals and Casimir force per unit area 
F (a) normalized to the force per unit area, 

acting between ideal metals are presented. 

As is seen in Fig. I, both the relative interaction energy and force between graphene and Au are greater than 
between graphene and Si at all separations. With the decrease of separation these quantities decrease. This is 
because E(a) and F(a) go to infinity more slowly than the respective dependencies for ideal metals in Eqs. (flU)) and 
(fTT)) . As a result, at the shortest separation indicated in Fig. I (a — 3nm) E/E = 0.0252, F/F — 0.0203 for Au and 
E/Eo = 0.0212, F/Fq = 0.0170 for Si. Note that at separations less than l-2nm the used model of graphene as a 
plasma sheet may become inapplicable because it does not take the atomic structure into account. With the increase 
of separation up to 1 /xm the interaction energy between graphene and Au achieves almost one half of that between 
ideal metals and the force magnitude is larger than 0.43 of that in the case of ideal metals. 

Now we compare the above numerical results for Au with the analytic calculations in the asymptotic region of short 
separations. The dielectric permittivity of Au along the imaginary frequency axis can be approximated by means of 
the plasma model, 

em = l + ^, (18) 

where ui p is the plasma frequency. At short separations only the TM mode in Eq. fTJ) contributes essentially to the 
van der Waals energy [2(|. Introducing the dimensionless variables v = aq and ( = a£/c and substituting Eq. (fT8|) in 
Eq. we rewrite the TM contribution to Eq. ([T]) in the following way: 

E(a) = E (a)f(a,p), (19) 



f(a, p) = T I 

~ Jo 



x In 



vdv J 



1 ^-^ A ^ e- 2v 



Here, Eq{o) was defined in Eq. (|16p and the two parameters are a = Sla and (3 — Lu p a/c. Instead of the variable £ we 
introduce the variable r\ = £/(/3v). This brings the correction factor / in Eq. (flU]) to the form 
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(20) 
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■ T 



l 



VP 2 



where i = a/ 1 . At short separations the parameters a and /3 are small (at the shortest separation a — 3nm it holds 
a = 2.02 x 10~ 3 , (3 = 0.137). It is easily seen that only small values of rj give the major contribution to the integral 
with respect to 77 in Eq. ([2U)l . Because of this, without loss of accuracy, one can replace the upper limit 1//3 for 00. 
Expanding the integrand in Eq. (|20p in powers of (3, one arrives at 



f(a,f3)=ph (t)+(3 3 h 1 (t) + 



where 



, , , 180 f°° 2 , f°° , 
h (t) = — / v dv drj 

I" JO JO 



(21) 



(22) 



-2v 



1 + 2f ) (1 + 2«V) 



x In 



Notice that the term of order /3 3 on the right-hand side of Eq. (|2Tj) and the respective terms originating from the 
TE contribution to Eq. ([1]) are small at short separations and can be neglected. The parameter t = a//3 2 is small 
and decreases with the increase of separation starting from the largest value of t = 0.108 at the shortest separation 
a = 3 nm. In Fig. 2 the values of ho are plotted as a function of t. 

As a result, the asymptotic representation of the van der Waals interaction energy between a graphene and an Au 
semispace at short separations is given by 



E(a) 



TP I \ UJ P a h 

c 



Qc 2 
u 2 a 



(23) 



The comparison of Eq. (j2"3"|) with the results of numerical computations in Fig. 1(a) (solid line) shows good agreement 
in the limits of 3% within the separation region from 15 to 80 nm. At the shortest separation of 3nm the error 
of the representation (|23[) achieves 11.5% if to compare with the computational results using the tabulated optical 
data. At the same time, within the separation region from 3 to 40 nm, Eq. (|23p is in agreement with the numerical 
computations using the plasma model (|18|) up to a maximal error of 3%. A few computational results are presented 
in Table 1. Column 1 contains separation distances, and in columns 2, 3 and 4 the values of the correction factor 
E(a)/ Eo(a) are contained computed by using the tabulated optical data, the plasma model dielectric function and 
the asymptotic representation (|23p . respectively. 

Note that the asymptotic expression (f2"3"]) at short separations contains the velocity of light, i.e., is relativistic. 
This is different from the van der Waals interaction at short separations of two bodies described by the dielectric 
permittivity [20] | but is in accordance with the short separation interaction of the two plasma sheets considered in 
Ref. 

It is instructive to consider the relationship between the graphene described by the model of a plasma sheet 
and the thin graphite plate containing a few hexagonal layers. For this purpose we compare the computed above 
interaction energy of graphene and Au semispace with the interaction of a graphite plate of some thickness d and 
the same semispace. The latter is computed using Eq. (fTJ) with the reflection coefficients © and ([5]). The dielectric 
permittivities of graphite along the imaginary frequency axis £ x ,y{i£,) and £ z («£) were computed in Ref. [18| by means 
of the Kramers-Kronig relations using the tabulated optical data for graphite from Ref. [34| . Here we use the results 
of that computations. 

In Fig. 3 the interaction energy density of the graphite plate and Au semispace normalized for the case of ideal 
metals in Eq. (fTB)) is plotted as a function of the relative plate thickness d/a at separations a = 1, 0.5, 0.3, 0.1, and 
0.05 jtzm (lines 1,2,3, 4, and 5, respectively). As is seen in Fig. 3, the magnitudes oiE/E decrease when d/a — > 0. This 
is in accordance with Eq. ([5]) for the reflection coefficients which go to zero when d vanishes. From this it follows that 
the graphene cannot be obtained from the macroscopic plate in the limit of zero thickness. There is, however, some 
relationship between the cases of graphene and the graphite plate. To trace it let us plot the values of graphene- Au 
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relative interaction energy at separations 1, 0.5, 0.3, 0.1, and 0.05 /im from Fig. 1(a) (the solid line) on the vertical 
axis in Fig. 3 and draw the tangents to the respective lines 1,2,3,4, and 5. It is seen that all tangents intersect 
with lines 1, . . . , 5 in the region d/a ~ 0.13 — 0.18. These values of the relative thickness are in fact characteristic 
for the graphite plates could be treated macroscopically in terms of the dielectric permittivity. As an example, at a 
separation of 50 nm (line 5 in Fig. 3) only the plates with thickness d > 0.15 • 50 = 7.5 nm are enough thick to be 
considered macroscopically. 



IV. INTERACTION OF GRAPHENE WITH A HYDROGEN ATOM OR A MOLECULE 



In Sec. II we have described the interaction of graphene with a material semispace by means of the Lifshitz-type 
formula ([T]) with the reflection coefficients Q and (fl~3| . This permits us to derive a Lifshitz-type formula for an atom 
near graphene by considering the limit of rarefied semispace [HI, [2(| • Let us expand the dielectric permittivity of a 
semispace in powers of the number of atoms per unit volume N preserving only the first order contribution [19J , 



e(i£) = 1 + 4na(i£)N + 0(7V 2 ), 



(24) 



where a(uj) is the dynamic polarizability of an atom. 

Expanding the reflection coefficients ([3]) and the energy density fl} in powers of N (see Ref. [26| for details), we 
arrive at the energy of the interaction between an atom and a graphene, 



E A (a) = 

V ; 2tt 



2r 



(i) 

TM,g 



k±dkj 

e 



d^a{i^)qe 



-2aq 



(25) 



q 2 c 2 



TE,g 



TM,j 



Recall that the reflection coefficients of graphene are given by Eqs. (fTB")) . 

As was shown in Ref. [Hj] , the atomic and molecular dynamic polarizabilities of H can be represented with sufficient 
precision using the single oscillator model, 



9a 



9m 



(26) 



Here, g a = a a (0)u> 2 is expressed through the static atomic polarizability a a (0) — 4.50 a. u. and the characteristic 

For a hydrogen molecule g m = a m (0)w 2 a where a m (0) = 5.439 a. u. and uj m = 14.09 eV 



frequency u a = 11.65 eV 

[HI . Note that before the substitution in Eq. ([25]) the polarizabilities (|26|) should be expressed in cubic meters taking 
into account that 1 a.u. of polarizability = 1.482 x 10~ 31 m 3 . 

For convenience in numerical computations, we introduce the dimensionless variable y = 2aq and represent Eq. 
in terms of the van der Waals coefficient C3: 



E A (a) 



C 3 (o) 



(27) 



h f°° 

= tk~ I d y e ~ 

16tt J 



cy/(2a) 



d£a(i£) 



2v 2 r (1) 

z y ' TM,g 



4a J t 





2c2 



,(1) 
TE.g 



TM,g 



The computational results for C3 as a function of a obtained by using Eqs. (j26|) and (p7| are presented in Fig. 4. The 
solid and dashed lines are related to the interaction of graphene with hydrogen atom and molecule, respectively. As 
is seen in Fig. 4, the magnitudes of the van der Waals coefficient for the hydrogen molecule interacting with graphene 
are larger than for the atom at all separations. Note that 1 a.u. of C3 = 0.646 x 10 -48 J m 3 . The comparison with the 
results of Ref. [HJ for the interaction of H atoms and molecules with graphite semispace shows that the interaction 
with graphene is by a factor of more than 2.5 weaker. 
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V. INTERACTION OF SINGLE- WALL CARBON NANOTUBES WITH METAL OR DIELECTRIC 

PLATE 



Let the single- wall nanotube of radius R lie along the y axis at a separation a from the boundary surface of a material 
semispace. For sufficiently small a the interaction energy in such configuration can be approximately obtained by 
using the proximity force theorem [36j and the Lifshitz-type formula fl} for the van der Waals energy between a 
graphene and a material semispace. According to the proximity force theorem, we replace the cylindrical surface by a 
set of infinitely long plane strips of width dx. The interaction between each strip, substituting a part of the cylindrical 
surface, and the opposite strip belonging to the boundary plane of a semispace is calculated by using Eq. JT]). The 
separation distance between the two opposite strips with coordinate x is 



z(x) 



R-y/R 2 



(28) 



Expanding the logarithms in Eq. Jl} in a power series, we present the interaction energy density between the strips 
in the form 



h f°° r 

\( r K 



1) (2) 
' TM r TM 



'-^ n 

n=l 



te'te 



(29) 



-2z(x)qn 



Recall that the reflection coefficients in the case under consideration are given by Eqs. and (flU)) . 

To find the interaction energy E n per unit length between the semispace and the nanotube, we integrate Eq. (|29p 
from x = — R to x = R (this is equal to twice the integral from zero to R). In so doing we replace the variable k± 
with q from Eq. ([2]): 



E n (a) 



h 

2^ 



(M) _(2) 
( TM TM 



qdq 



cq oo 

d£ V -e- 2aqn 



r (l) r (2) 
'TE'TE 



n=l 

■R. 



(30) 



dxe 



-2qn(R-VW 



By introducing the new variable s = 1 — yl — x 2 /R 2 , the integral / with respect to x in Eq. ([30]) can be written in 
the form 



I = R 



1 (I - s)ds 
o V s (2-s)' 



~2qnRs 



(31) 



The major contribution in Eqs. (|30[) and pip comes from q ~ 1/a. Bearing in mind that the proximity force theorem 
works good for R 3> a, we conclude that the magnitude of the integral / is determined by the behavior of the integrand 
around the lower integration limit. Neglecting s compared with unity in Eq. (|31[) . we arrive at 



E> /■! e ~2qnRs 




(32) 



where Erf(z) is the error function. Using once more the conditions q ~ 1/a and R 3> a, we conclude that Erf(^/2gni?) 
1 and obtain from Eq. (|30"|) 



E n (a) 



At: 2 



7o Jo £=J n^/n 
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r I TE TE 



(33) 
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'Air 3 / 2 



j Vqdq / dt, 
Jo Jo 



Li 



3/2 ( r TM r l 



2) „-2a 9 
TM e 



Li 
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(r TE r. 
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TE e 



where Li 3/ / 2 (z) is the polylogarithm function. 
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By analogy, for the force per unit length it follows 



|^ r^dqTdi (34) 



X 



Li 



1/2 \^ TM TM J + L1 l/2 ^ r TE r TE e J 



Note that for an ideal metal cylinder in close proximity to an ideal metal plate the interaction energy and force per 
unit length are given by [37| 

E5(*) = -\\f2^ t py(a)=-lJ?J^. (35) 
a 2 V a 960V2 a 3 V a 384^2 

It is significant that Eqs. ([55)) and (f3"4"| follow also in the limit of short separations from the exact result for a 
nanotube described by the cylindrical surface with the boundary conditions (fTU)) located above a material semispace. 
This result can be obtained by the method of functional determinants used previously in Refs. [13, [H, Hi| for the 
configurations of a sphere or a cylinder made of ideal metal placed above an ideal metal plate (the details will be 
published elsewhere). The exact result contains also the asymptotic case of large separations n» Ji between the 
nanotube and the metal semispace described by the plasma model (1181) : 



he r°° \/(— ) +P 2 ~P 



W ' J ° y( t ) +p 2 +p 

This asymptotic behavior does not depend on the nanotube parameter J7. In the limit uj p — * 00 the integral in Eq. (|36[) 
tends to 1 /4 and we arrive at the previously known result 

£o» = s % b. ' (37) 
B7rcr in ^ 

valid for an ideal metal cylinder far away from an ideal metal plate (37l . [38| ■ 

Now we present the computational results for a nanotube in close proximity and far away from a semispace. 
Consider, first, the nanotube of radius R in close proximity to Au or Si semispace and compute the normalized 
interaction energy E n /Eq and the force F n /F n using Eqs. (f33|) (f35|) . As in Sec. Ill, the dielectric permittivities of 
Au and Si along the imaginary frequency axis are found from the Kramcrs-Kronig relations and tabulated optical 
data for the complex index of refraction. A few computed results are presented in Table II. Column 1 contains 
the separation distance, columns 2 and 3 contain the values of the correction factor E n (a)/EQ (a) for Au and Si, 
respectively, and columns 4 and 5 contain the analogous correction factor to the force F n (a) / Fy (a) . As is seen in 
Table II, the magnitudes of all correction factors for Si are smaller than for Au and they are monotonously increasing 
functions with the increase of separation. Note that the normalized values in Table II do not depend on the nanotube 
radius. However, bearing in mind that the largest diameter of a single- wall carbon nanotube is of about 10 nm Q], 
the separation region where the approximate equations (|33p and (|34[) arc applicable is very narrow. 

Now we consider the nanotube of radius R far away from an Au surface described by the plasma model (b» R). In 
this case the Casimir interaction energy is given by Eq. (I36p . The results of numerical computations of the normalized 
interaction energy E n (a)/EQ(a), where the energy Eq(o) for an ideal metal case is given in Eq. (|3T|) . are presented in 
Fig. 5 as a function of separation. As is seen in Fig. 5, the normalized interaction energy increases with the increase 
of separation approaching the case of an ideal metal cylinder above an ideal metal plate. At a = 1 /im it achieves 
the value of 0.956 very close to the ideal metal case. For the typical single-wall nanotube diameter of about I nm the 
asymptotic expression (|36p is applicable in the wide separation region from a few nanometers to about I micrometer 
(for larger separations the thermal corrections should be taken into account). 



VI. CONCLUSIONS AND DISCUSSION 



In the above we have obtained Lifshitz-type formulas describing the van der Waals and Casimir interaction between 
graphene and material plate, between an atom or a molecule and graphene, and between a single- wall carbon nanotube 
and material plate. The distinguishing feature of these formulas is that they describe graphene by using the reflection 
coefficients obtained from the specific boundary conditions for the electromagnetic oscillations on the infinitely thin 
plasma sheet. This permits to circumvent the use of the concept of dielectric permittivity commonly used in the Lifshitz 
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theory of the van der Waals and Casimir force between macrobodies, but being not directly applicable to single-wall 
carbon nanostructures. The developed formalism is supplementary to widely applied theoretical approaches describing 
the interactions of carbon nanostructures with macrobodies, atoms or molecules, e.g., to the density functional theory 
(see Introduction). Together with Ref. [l|j], where the Lifshitz theory was applied to the multi-wall carbon nanotubes, 
it provides the foundation for the application of quantum statistical physics to the investigation of dispersion forces 
in carbon nanostructures. 

The obtained Lifshitz-type formulas for the van der Waals and Casimir energy and force were applied to the case of 
graphene interacting with Au and Si walls. The wall material was described by the dielectric permittivity along the 
imaginary frequency axis computed using the tabulated optical data for the complex index of refraction for Au and 
Si. In the case of Au the analytic asymptotic expression for the interaction energy with graphene was also obtained 
using the plasma model dielectric function. The relationship between the graphene, described by the plasma sheet, 
and the thin graphite plate was investigated. As an example of microparticle-carbon nanostructure interaction, we 
have calculated the van der Waals coefficients for the interaction of hydrogen atoms and molecules with graphene. 
Another example, considered in the paper, is the van der Waals or Casimir interaction of a material wall with a 
single-wall carbon nanotubc in close proximity or far away from the wall, respectively. As was noted in Sec. Ill, the 
obtained results for the interaction energy and force are applicable for practical calculations at separations larger than 
1-2 nm. At smaller separations there may be the attractive forces of chemical nature. At separations less than 1 nm 
short-range repulsive forces of exchange nature come into play. These forces depend on atomic structure of a surface 
and cannot be described macroscopically by means of the boundary conditions. At intermediate separations between 
the exchange repulsion and van der Waals attraction some phenomcnological potentials can be used for practical 



In the future work it would be interesting to perform the comparative computations of the van der Waals interaction 
with carbon nanostructures by using different theoretical approaches. The suggested application of the Lifshitz theory 
to graphene and single-wall carbon nanotubes can be extended to more complicated structures such as fullerene 
molecules or graphitic cones The interaction of hydrogen atoms and molecules with the single- wall carbon 

nanotubes can be also considered on the same footings. This can be done approximately by using the proximity force 
theorem because the immediate application of the Lifshitz-type formulas is possible for only planar structures (the 
results of this work in progress will be published elsewhere). In this paper we did not deal with the thermal effects. 
However, the generalization of the proposed formalism to the case of nonzero temperatures is straightforward. 
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Figures 




FIG. 1: The normalized to the case of ideal metals van der Waals and Casimir energy (a) and force (b) per unit area between 
a graphene and a semispace versus separation. The solid and dashed lines are related to the semispace made of Au and Si, 
respectively. 
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FIG. 2: The function ho in Eq. (|21|l versus the dimensionless parameter t = Q.c? /(ajpO). 



15 




FIG. 3: The normalized to the case of ideal metals van der Waals and Casimir energy per unit area between a graphite plate 
of thickness d and Au semispace versus relative thickness. Lines 1,2,3,4, and 5 are related to separations a = 1, 0.5, 0.3, 0.1, 
and 0.05 fj,m, respectively. The respective normalized interaction energies of a graphene with an Au semispace are marked on 
the vertical axis (see text for further discussion). 
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FIG. 4: The van der Waals coefficient C3 for the interaction of a hydrogen atom (the solid line) and a molecule (the dashed 
line) with graphene versus separation. 
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FIG. 5: The normalized to the case of ideal metals Casimir energy per unit length between a carbon nanotube and Au semispacc 
versus separation. 
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Tables 
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TABLE I: Ratios E(a)/ Eo(a) for the Casimir interaction of an Au semispace and graphene computed using: (a) — tabulated 
optical data for Au, (b) — plasma model for Au, (c) — asymptotic expression (|23[1 , 



a 




E(a)/E (a) 




M 


(a) 


GO 


(c) 


3 


0.0252 


0.0222 


0.0223 


5 


0.0338 


0.0307 


0.0310 


10 


0.0496 


0.0466 


0.0473 


15 


0.0618 


0.0590 


0.0602 


20 


0.0720 


0.0694 


0.0711 


30 


0.0893 


0.0870 


0.0895 


40 


0.104 


0.102 


0.105 


50 


0.117 


0.115 


0.119 


60 


0.129 


0.127 


0.131 


70 


0.139 


0.138 


0.143 


80 


0.1496 


0.1485 


0.154 



20 



TABLE II: Ratios E n (a)/EQ(a) and F n (a)/FJ l (a) for the van der Waals interaction of a carbon nanotube with Au and Si 
semispaces 



a 

(nm) 


Au 


E n (a)/ES(a) 


o • 
Bl 


Au 


F»(o)/jy(o) 


o • 
Bl 


1 


0.0151 




0.0126 


0.0114 




0.00945 


1.5 


0.0193 




0.0162 


0.0147 




0.0123 


2 


0.0230 




0.0193 


0.0175 




0.0147 


2.5 


0.0262 




0.0221 


0.0201 




0.0169 


3 


0.0291 




0.0245 


0.0224 




0.0189 



